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The strength functions of quadrupole modes in the unstable oxygen isotopes 22 and 24 are 
calculated using an extended version of the time-dependent Hartree-Fock theory known as the time- 
dependent density-matrix theory (TDDM). It is found that TDDM gives the lowest quadrupole 
states which are energetically shifted upward and become significantly collective due to the coupling 
to two-body configurations. It is pointed out that these features of the lowest quadrupole states are 
similar to those obtained in the quasi-particle random phase approximation. 

PACS numbers: 21.10.Re, 21.60.Jz, 27.30.+t 




The study of low-lying quadrupole states in unstable 
oxygen isotopes has recently gained both experimental 
§, §, |, f§ and theoretical interest §, §, @, §, §. One 
of the reasons for this is that the energy and transi- 
tion strength of the first 2 + state are considered to be 
closely related to the change in the shell structure in this 
region of neutron rich nuclei |j~T[] . Various theoret- 
ical approaches have been used to study the low-lying 
quadrup_ole_states in unstable oxygen isotopes: The shell 
the random-phase approximation (RPA) 
and the quasi-particle RPA (QRPA) @ |, | . 
It is shown by the QRPA calculations |t], |f, |jj that paring 
correlations are essential to reproduce the large B(E2) 
values of the first 2 + states. Earlier TDDM calculations 
p5[ did not include the spin-orbit splitting of the single- 
particle states, thereby the first 2 + states, which mainly 
consist of inner shell transitions of neutrons, could not be 
treated properly. The TDDM approach has recently been 
improved to include the spin-orbit force and now is able 
to handle low-lying quadrupole states in a more quanti- 
tative way [[l6]]. In this letter we apply TDDM to study 
low-lying quadrupole states, especially, the first 2 + states 
and show that TDDM gives results similar to QRPA. The 
present formulation of TDDM assumes the Hartree-Fock 
(HF) ground state as a starting ground state as will be 
shown below. The quantitative application of TDDM, 
therefore, is restricted to the quadrupole states in 22 O 
and 2 O, where the first-order approximation of the sub- 
shell closure may be justified |12j . 

TDDM is an extended version of the time-dependent 
Hartree-Fock theory (TDHF) and is formulated to deter- 
mine the time evolution of one-body and two-body den- 
sity matrices p and p2 in a self-consistent manner [fj"7| . 
The equations of motion for p and p2 can be derived 
by truncating the well-known Bogoliubov-Born-Green- 
Kirkwood-Yvon hierarchy for reduced density matrices 
|fl8| . To solve the equations of motion for p and p2, p 
and C2 (the correlated part of p^) are expanded using a 
finite number of single-particle states ip a which satisfy a 



TDHF-like equation, 

p(U',t) = £n QQ <(i)</>a(l,t)C'(l'.*), (1) 

C 2 (121'2',*) = P2-A{pp) 

= ^ C a f3 a >p>(t) 

x Wpfr W a > (!', ( 2 '. *). ( 2 ) 

where A is the antisymmetrization operator and the num- 
bers denote space, spin, and isospin coordinates. Thus, 
the equations of motion of TDDM consist of the following 
three coupled equations M|: 



ih^-Jj a (l,t) = h(l,t)^ a (l,t), 
ot 



(3) 



ihh aa > = '^2[(a(3\v\"fS)C 7 sa'ii - C al3l shS\v\a' f3)}, (4) 



ihCapa'p' — B a p a 'f3> + P a 0a'/3' + Hafia'P 1 



(5) 



where h is the mean-field Hamiltonian and v the residual 
interaction. The term B a fia'$' on the right-hand side of 
Eq.(5) represents the Born terms (the first-order terms 
of v). The terms P a p a 'f3' and H a a ipi in Eq.(5) con- 
tain G a p a ipi and represent higher-order particle-particle 
(and hole- hole) and particle-hole type correlations, re- 
spectively. Thus full two-body correlations including 
those induced by the Pauli exclusion principle are taken 
into account in the equation of motion for C a p a i pi . The 
explicit expressions for B a p a ip>, P a p a 'p' an d Hafta'p' are 
given in Ref . Jl7t . Diverse configurations including ex- 
citations of the 16 O core can be taken into account in 
Eqs.(3)-(5). In order to obtain a clear understanding of 
contributions of various configurations we solve the cou- 
pled equations in the following three schemes: 1) Two- 
body amplitudes C a p a ip> consisting of the neutron 2s 1/2, 
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ld 3 /2 and ld 5 / 2 states only are considered. We call this 
scheme TDDM1. 2) The neutron 2p 3 / 2 an d I/7/2 are 
added to TDDM1. We call this version TDDM2. 3) 
All single-particle states up to the 2p 3 / 2 and I/7/2 are 
taken into account for both protons and neutrons. We 
call this scheme TDDM3. Two-body configurations cor- 
responding to excitations of the 16 O core are included in 
TDDM3. The number of independent C a (3 Q '^'s drasti- 
cally increases from 513 (TDDM1) to 50,727 (TDDM3). 

The E2 strength function is calculated according to 
the following three steps : 

1) A static HF calculation is performed to obtain the 
initial ground state. The ld 5 / 2 an d % s i/2 orbits are as- 
sumed to be the last fully occupied neutron orbit of 22 O 
and 24 0, respectively. The Skyrme III with spin-orbit 
force is used as the effective interaction. The Skyrme III 
has been used as one of standard parameterizations of 
the Skyrme force in nuclear structure calculations even 
for very neutron rich nuclei |t], The single-particle 
wavefunctions are confined to a cylinder with length 20 
fm and radius 10 fm. (Axial symmetry is imposed to cal- 
culate the single-particle wavefunctions p0||.) The mesh 
size used is 0.5 fm. 

2) To obtain a correlated ground state, we evolve the 
HF ground state using the TDDM equations and the fol- 
lowing time-dependent residual interaction 



v(t) = (l-e-*/ T )i;(r-r'). 



(0) 



The time constant r should be sufficiently large to ob- 
tain a nearly stationary solution of the TDDM equations 
plj . We choose r to be 150 fm/c. In a consistent calcu- 
lation the residual interaction should be the same as that 
used to generate the mean field. However, a Skyrme-type 
force contains momentum dependent terms , which make 
the computation time of two-body matrix elements quite 
large. Therefore, we need to use a simple force of the S 
function form v oc S 3 (r ~ r'). In order to make a com- 
parison with the results of HFB and QRPA calculations, 
we use the following pairing-type residual interaction of 
the density-dependent 5 function form |23] 



v(r - r') = u (l - p(r)/ Po )S 3 (r - r'), 



(7) 



where p{r) is the nuclear density, po and vq are set to 
be 0.16fm~ 3 and —1200 MeV fm 3 , respectively. Similar 
values of po and vq have been used in the Hartree-Fock- 
Bogoliubov (HFB) calculations p3, [2|, ||| in truncated 
single-particle space. The time step size used to solve the 
TDDM equations is 0.75 fm/c. 

3) The quadrupole mode is excited by boosting 
the single-particle wavefunctions at t = 5r with the 
quadrupole velocity field: 



V>a(5r) 



JkQ 



1p a . 



(8) 



We consider the isoscalar quadrupole mode and the pro- 
ton quadrupole mode to calculate a B(E2) value. For 
the former, Q is equal to r 2 Y 2 Q for both protons and neu- 
trons, and for the latter only the proton single-particle 



states are boosted. When the boosting parameter k is 
sufficiently small, the strength function defined by 



S(E) ^J2\^\Q^o)\ 2 5(E - E n 



(9) 



is obtained from the Fourier transformation of the time- 
dependent quadrupole moment Q(t) as 



S(E) 



1 

irkh 



Et 

Q(t) sin -^-dt. 



(10) 



In Eq.(9) |$o) is the total ground-state wavefunction and 
is the wavefunction for an excited state with exci- 
tation energy E n , and Q the quadrupole operator. The 
TDDM calculations are stopped at t — 1500 fm/c and the 
upper limit of the time integration in Eq.(10) becomes 
750 fm/c. To reduce fluctuations in S(E), originating 
from the finite-time integration, the quadrupole moment 
is multiplied by a damping factor e r vt l 2h with T — lMeV 
before the time integration. Since the integration time is 
limited, the strength function in a very low energy region 
(E < 2nh/750 « 2MeV), which depends on long time 
behavior of Q(t), is not well determined and somewhat 
affected by numerical inaccuracies. The energy-weighted 
sum rule (EWSR) for the isoscalar quadrupole mode is 
expressed as 



S(E)EdE 



±(%\[Q,[H,Q]]\<S> ) 



5ft 2 

47T77J 



A(r 2 



(11) 



where A is the mass number, m is the nucleon mass and 
(r 2 ) means the mean square radius of the ground state. 
For the proton mode the mass number is replaced by the 
atomic number and the mean square radius is taken for 
protons, and there appears an enhancement factor due 
to the momentum dependence of the Skyrme interaction 
which is expressed as 



r 2 p p (r)p n (r)dr, 



(12) 



where t\ and t 2 are the parameters of the momentum 
dependent parts of the Skyrme force, and p p (p n ) is 
the proton (neutron) density distribution. In TDDM3 
where proton single-particle states are also considered 
in two-body configurations, a term containing C a j3 a ipi, 
which exhibits the effects of ground-state correlations 
appears in addition to Eq.(12) jig . However, its con- 
tribution is less than 1% of the total EWSR value and 
is neglected. The contribution of Eq.(12) is about 15% 
of the total EWSR value. The strength function ob- 
tained in TDHF (without Eqs.(4) and (5)) is equivalent 
to that in RPA without any truncation of unoccupied 
single-particles states because the TDHF equations for 
the boosted single-particle wavefunctions ip a is solved in 
coordinate space. The boundary condition for the con- 
tinuum states, however, is not properly taken into ac- 
count in our calculation because all the single-particle 
wave functions are confined to the cylindrical geometry. 
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We first discuss some ground state properties of O 
and 24 evaluated at t = 5r in TDDM1 where only the 
neutron 1^5/2, 2si/2 and 1^3/2 orbits are use to evalu- 
ate C a p a ipi. In 22 O, the energies of the neutron 1 0(5/2, 
2si/ 2 , and ld 3 / 2 states are — 6.8MeV, — 2.3MeV, and 
1.3MeV, respectively. Their occupation probabilities are 
0.92, 0.14, and 0.05, respectively. In 24 0, their ener- 
gies and occupation numbers are — 7.3MeV, — 2.3MeV, 
0.8MeV, 0.97, 0.89, and 0.10, respectively. The correla- 
tion energy E COI defined by 

E ™ = ^<a/?M 7 <5)C 7fo/3 (13) 

is -4.4 McV in 22 and -3.5 MeV in 24 0. The correla- 
tion energy defined above includes effects of various two- 
body correlations. However, pairing correlations seem to 
dominate. This is because E COT does not change much 
even if the sum over single-particle states in Eq.(13) is 
restricted to a pair of conjugate states under time rever- 
sal such as 

iX>/3M7<5)C 7fa/3 — > i£<aaM/?j8)Cag aaj (14) 

where a is a state conjugate to a under time reversal. 
The correlation energies thus obtained for 22 O and 24 O 
are —4.8 MeV and —3.9 MeV, respectively. Since the 
pairing correlations seem dominant, it might be meaning- 
ful to estimate a quantit y A corresponding to an average 
pairing gap in HFB p^, pjfl. We define A as 

where v a and u a are calculated as y/n aa and \/T— n aa , 
respectively. The values of A thus obtained are —3.1 
MeV and -2.7 MeV for 22 O and 24 O, respectively. These 
values may be comparable to the empirical pairing gap of 
A « 12/a/A w 2.5 MeV. When the single-particle space 
is extended to the 2p 3/2 and l/ 7/2 in TDDM2, A be- 
comes slightly larger and its difference between 22 O and 
24 becomes smaller: The values of A in TDDM2 are 
-3.6 MeV and -3.5 MeV for 22 and 24 0, respectively, 
and become close to the values obtained from HFB calcu- 
lations using a Woods-Saxon potential ||. These values 
might be rather large as compared to experimental gaps 
obtained by combining the binding energies of neighbour- 
ing nuclei. However, only a semi-quantitative agreement 
with experimental values should be expected because A 
in Eq.(15) has no direct relation to the definition of the 
experimental gap. 

The strength functions for the isoscalar quadrupole 
modes in 22 O and 24 O are shown in Figs.l and 2, respec- 
tively. The dotted lines depict the results in TDDM1 and 
the solid lines those in TDDM2. The results in TDHF 
are also shown with the dot-dashed lines. The width of 
each peak is due both to the smoothing with r = 1 MeV 
and to the finite-time integration of Eq.(10). In the very 
low energy region (E < 2 MeV), S(E) fluctuates slightly 
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FIG. 1: Strength distributions of the isoscalar quadrupole 
modes in 22 O calculated in TDHF (dot-dashed line), TDDM1 
(dotted line) and TDDM2 (solid line). 
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FIG. 2: Strength distributions of the isoscalar quadrupole 
modes in 24 calculated in TDHF (dot-dashed line), TDDM1 
(dotted line) and TDDM2 (solid line). 

and has small negative components. This originates in 
uncertainties associated with the calculation of Eq.(10) 
as mentioned above. The fraction of the EWSR values 
depleted below 60MeV is about 98% in all calculations 
shown in Figs.l and 2. Based on the single-particle en- 
ergies of the neutron 2s and Id states, we consider that 
the lowest state seen at E = 3.8 MeV in TDHF for 22 O 
originates in the transition of a neutron from the ld 5 / 2 
orbit to the 2si/2 orbit, while that at E = 2.3 MeV in 
24 comes from transition of a neutron from the 2s!/ 2 
orbit to the 1^3/2 orbit. The second lowest states seen 
at E = 7.8 MeV in the TDHF calculations for both nu- 
clei are also considered to be from the inner shell transi- 
tion of a neutron, that is, the transition from the ld§/2 
to the ld 3 / 2 - When the two-body correlations among 
the neutron ld 5 / 2 , 2si/ 2 and ld 3 / 2 states are included in 
TDDM1, the collectivity and excitation energies of the 
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FIG. 3: Strength distributions of the proton quadrupole 
modes in 22 O calculated in TDHF (dot-dashed line), TDDM1 
(dotted line) and TDDM2 (solid line). 

lowest states are increased and the strength of the gi- 
ant quadrupole resonance (GQR) located at 20 MeV is 
slightly decreased correspondingly. The increase of the 
excitation energy of the lowest state is from 3.8 MeV to 
4.7 MeV in 22 and from 2.3 MeV to 3.5 MeV in 24 0. 
The increase of the collectivity of the lowest states is 
due to the mixing of quadrupole states consisting of two- 
body configurations and the upward shift in excitation 
energy is related to the lowering of the ground state due 
to two-body correlations. The paring correlations seem 
to dominate the two-body correlations because the two- 
body amplitudes of C a& as type are most important as 
discussed above. These behaviors of the first 2 + states 
under the influence of the pairing correlations are sim- 
ilar to those obtained in QRPA @, §. In TDDM the 
pairing correlations are treated as two-body correlations, 
whereas they are considered in a mean field approxima- 
tion in QRPA. The expansion of the single-particle space 
up to the neutron 2p 3 / 2 and I/7/2 states (TDDM2) fur- 
ther enhances the collectivity and the excitation energy: 
The excitation energy of the lowest state is increased to 
5.0 MeV in 22 and to 4.4 MeV in 24 0. The quadrupole 
strengths below 6 MeV in 22 are 214 fm 4 , 271 fm 4 
and 300 fm 4 in TDHF, TDDM1 and TDDM2, respec- 
tively. The quadrupole strengths below 6 MeV in 24 O 
are 163 fm 4 , 206 fm 4 and 249 fm 4 in TDHF, TDDM1 
and TDDM2, respectively. The increase in the collectiv- 
ity with increasing number of single-particle states seems 
to be slightly larger in 24 than in 22 0. 

The proton quadrupole modes in 22 and 24 are 
shown in Figs. 3 and 4, respectively. The fraction of 
the EWSR values depleted below 60MeV is about 93% 
in all calculations shown in Figs. 3 and 4. The strength 
distributions seen above 25 MeV correspond to isovec- 
tor quadrupole resonances. The B(E2) values obtained 
from the integration of S(E) below 6 MeV in 22 O are 11.0 
e 2 fm 4 , 14.2 e 2 fm 4 and 16.4 e 2 fm 4 in TDHF, TDDM1 and 
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FIG. 4: Strength distributions of the proton quadrupole 
modes in 24 calculated in TDHF (dot-dashed line), TDDM1 
(dotted line) and TDDM2 (solid line). 

TDDM2, respectively. The B(E2) value of 16.4 e 2 fm 4 in 
TDDM2 may be comparable to the observed value of 
21 ± 8 e 2 fm 4 || for the first 2+ state, although the ex- 
citation energy of 5.0 MeV is larger than the experimen- 
tal value of 3.2 MeV @. The QRPA calculations §, § 
also give the lowest quadrupole states in oxygen isotopes 
whose excitation energies are slightly larger than exper- 
imental values. The excitation energies of the first 2 + 
states sensitively depend on the single-particle energies 
of the 2s and Id states. To reproduce the observed ex- 
citation energies, appropriate choice and adjustment of 
the Skyrme force parameters including the strength of 
spin-orbit force may be necessary S. It has also been 
pointed out |2(| that the energy of the 2si/ 2 state in oxy- 
gen isotopes is significantly shifted downward due to the 
coupling to phonon states. The B(E2) values obtained 
from the integration of S(E) below 6 MeV in 24 are 3.5 
e 2 fm 4 , 5.4 e 2 fm 4 and 7.2 e 2 fm 4 in TDHF, TDDM1 and 
TDDM2, respectively. The small B{E2) value in 24 is 
variant from the QRPA calculation done by Matsuo || 
but consistent with other QRPA calculations ^ |j| and 
the shell model calculations [|[ |). 

The strength function of the isoscalar quadrupole 
modes in 24 O calculated in TDDM3 (solid line) where 
the excitations of the 16 O core are included is shown in 
Fig. 5. To obtain a stable ground state in TDDM3, we 
reduced the strength of the residual interaction and used 
v = -800 fm 3 MeV in stead of v = -1200 fm 3 MeV. 
The result in TDDM2 calculated with this value of vq 
is also shown in Fig. 5 with the dotted line. The highest 
peak of GQR is reduced and GQR becomes broader in 
TDDM3. Thus the major effect of the excitations of the 
16 O core on the quadrupole mode is to modify the dis- 
tribution of the GQR strength. However, the damping 
of GQR in 24 is modest as compared with that in 16 0. 
To make this point clear, we show in Fig. 6 the strength 
function for the isoscalar quadrupole mode in 16 O. The 
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solid line denotes the result in TDDM3 with vq = —800 
fm 3 MeV and the dotted line that in TDHF. The EWSR 
value in TDDM3 is 7 % larger than that in TDHF. The 
splitting of the GQR strength obtained in TDDM3 is con- 
sistent with experimental observation [ p7| . The damping 
of GQR in 16 O is caused by its coupling to two-body 
configurations mainly consisting of two holes in the lp 
states and two particles in the 2s or Id states. In 24 0, 
the coupling of GQR to two-body configurations is hin- 
dered by the spatial extension of neutron single-particle 
wavefunctions due to the neutron excess. 



FIG. 5: Strength distributions of the isoscalar quadrupole 
modes in 24 O calculated in TDDM3 (solid line) and TDDM2 
(dotted line). The strength of the residual interaction is no = 
-800 fm s MeV. 
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FIG. 6: Strength distributions of the isoscalar quadrupole 
modes in 16 O calculated in TDDM3 (solid line) and TDHF 
(dotted line) . The strength of the residual interaction used in 
TDDM3 is v = -800 fm 3 MeV. 



In summary, the strength functions of the quadrupole 
resonances in 22 and 24 were studied using TDDM. In 
this approach, the correlated ground state was first ob- 
tained starting from the HF ground state with the sub- 
sequent excitation of the quadrupole mode. It is shown 
that the lowest quadrupole states are shifted upward and 
become significantly collective due to the coupling to 
two-body configurations. The obtained B(E2) value for 
22 O was found comparable to the observed one. It was 
pointed out that the lowest quadrupole states obtained in 
TDDM have properties similar to those in QRPA. It was 
also found that the damping of GQR in 24 is modest as 
compared with that in 16 O. 
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